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The  scattering  of  plane  sound  waves  from  an  air-filled  steel  spherical  shell  submerged  in  water 
in  the  frequency  band  0<A:,a<500  is  studied.  This  analysis  is  based  on  a  methodology  [Ayres 
etal.,  Int.  J.  Solids  Struct.  23, 937-946  (1987)  and  G.  Gaunaurd  and  M.  F.  Werby,  J.  Acoust. 

Soc.  Am.  82,  2021-2033  ( 1987)  ]  proposed  that  uses  the  exact  three-dimensional  equations  of 
dynamic  elasticity  to  describe  the  shell  motions  and  to  predict  its  sonar  scattering  cross  section. 

This  approach  is  valid  at  all  frequencies,  for  shells  of  any  thickness,  of  any  (constant) 
curvature,  and  it  accounts  for  their  fluid-loaded  condition.  The  methodology  is  used  to  predict 
the  cross  sections,  which  are  later  interpreted  on  the  basis  of  the  various  resonance  features 
that  manifest  themselves  in  the  frequency  response.  The  spectral  locations  of  these  resonances 
depend  on  the  various  types  of  elastic  waves  propagating  along  the  shell,  or  in  the  surrounding 
fluid.  The  exact  plots  are  genet  ated  for  the  phase  {cp )  velocities  of  these  (Lamb)  waves  always 
accounting  for  the  curvature  and  fluid-loading  effects  present  on  the  shell,  without  appeals  to 
plate  waves  or  theories.  Some  of  the  dispersion  plots  were  generated  using  the  Donnell  shell- 
theory  approximation,  which  seems  tc  yield  accurate  results  up  to  the  coincidence  frequency. 

Aside  from  the  broad  resonance  lobe  present  at  the  coincidence  frequency,  there  is  anotht  e 
high-frequency  resonance  lobe,  due  to  a  thickness-resonance  effect,  which  was  also  predicted 
and  displayed.  A  partial-wave  analysis  of  the  resonance  response  curve  for  a  thin  shell,  around 
its  coincidence  frequency,  serves  to  identify  the  origins  of  the  various  types  of  observed 
resonance  features  and  to  relate  them  to  the  elastic  and  acoustic  waves  that  propagate  along 
the  shell  or  the  outer  fluid.  Many  computer-generated  graphs  are  displayed  to  illustrate  the 
above  points. 
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INTRODUCTION 

The  study  of  the  reflection  and  scattering  of  sound 
waves  from  submerged  elastic  shells  «.  .nstitutes  a  problem 
area  not  only  very  rich  in  scientific  and  technical  challenges 
but  also  of  crucial  importance  to  our  employer  in  view  of  its 
connection  with  the  sensing  and  the  identification  of  sub¬ 
merged  scatterers.  Work  on  this  subject  started  proliferating 
in  the  early  1960s,1'3  and  over  the  years  it  has  included 
bare1-3  and  coated3  shells  of  spherical1"3  and  cylindrical4 
geometries,  among  others. 

The  methodology  that  we  will  use  here  to  generate  our 
predictions  and  interpretations  is  a  three-dimensional  elasti¬ 
city  approach  that  first  appeared5  6  in  1987.  We  will  consid¬ 
er  Ref.  6  as  the  first  part  of  this  work,  hence,  its  almost 
identical  title.  We  have  worked  on  this  subject  over  the 
years,7"10  using  either  classical  approaches  or  specific  reso¬ 
nance  techniques  particularly  useful  in  the  (broad)  reso¬ 
nance  region  of  submerged  elastic  structures.  We  have  sum¬ 
marized  many  of  our  findings*"10  in  reviews  that  extensively 
cite  the  truly  large  number  of  contributions  to  this  important 
area,  from  all  over  the  world.  Of  particular  importance  are  a 
series  of  monographs11"16  and  papers17,13  that  deal  with 
Rayleigh  and  Lamb  waves.  It  has  become  undisputed  over 


the  years  that  as  sound  waves  impinge  on  elastic  scatterers, 
the  returned  echoes  have  features  that  depend  on  the  type  of 
elastic  waves  that  get  excited  on  the  scatterer.  For  either  an 
elastic  half-space  or  for  an  elastic  layer  of  finite  thickness,  the 
basic  waves  are,  respectively,  the  Rayleigh  and  the  (various 
types  of)  Lamb11  waves.  For  convex  solid  (or  hollow)  elas¬ 
tic  t  odies  immersed  in  fluids  these  waves  have  retained  their 
names  >or  the  flat  case,  but  have  been  generalized1317  to 
account  for  the  nonvanishing  body  curvature.  It  is  in  this 
context  that  these  names  are  currently  used.  In  spite  of  the 
fact  that  exact  three-dimensional  descriptors  of  the  shell  mo¬ 
tions  are  available,3,5,6  the  tendency  to  explain  features  in  the 
sonar  cross  sections  of  submerged  shells  in  terms  of  (loaded 
or  unloaded )  flat  plate  waves  or  theories,  still  seems  to  con¬ 
tinue.  We  use  the  predictions  of  our  approach;5,6  (i)  to  ex¬ 
plain  a  variety  of  new  features  that  have  emerged,  and  (ii)  to 
compare  them  in  some  instances  to  the  results  obtained  from 
the  above  mentioned  simpler  models. 

I.  THEORETICAL  BACKGROUND 

The  complete  classical  and  RST  formulation  required 
for  the  study  of  the  acoustic  echoes  scattered  from  air-filled 
elastic  shells  in  water  has  been  given  in  Refs.  5  and  6  for  bare 
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spherical  sheLls,  in  Ref.  3  for  viscoelastically  coated  spherical 
shells,  in  Refs.  4  and  1 8  for  bare  cylindrical  shells,  and  in  Ref. 
19  for  viscoelastically  coated  cylindrical  shells.  The  essen¬ 
tials  of  this  formulation  arc  contained  in  the  form  function, 
which  is  given  by  the  expression 


I  /„  (*.*)!  = 


n  *-jO 


—  f  (2«+  1  )A„Pn  (cos  8) 
ix  n~0 


(1) 


where  x  =  kxa,  and  the  coefficients  A„  for  the  bare,  air- 
filled,  spherical  sh  ill  in  water,  are  given  by  the  ratios  of  two 
6x6  determinants  B„  and  D„,  viz.. 


A„(x)=Bn(x)/D„(x),  (2) 

where  the  elements  of  these  determinants  are  given  in  Refs.  5 
and  6.  The  outer  (or  inner)  radius  of  the  shell  is  a  (or  b). 

There  are  many  ways  to  proceed  from  here  on.  Some  of 
these  are  enumerated  below, 

(i)  Gene  ration  of  the  backscattering  or  monostotic  form 
function:  We  compute  and/or  display:  \fx  (8  =  tt,x)  |  vs  x. 
Here  we  use  P„  (cos  •<.-)  =  (  —  1 )". 

(ii)  Isolation  of  the  resonances — spectrogram:  Here  we 
compute  and/or  display  |/w  (6  =  tt,x)  —  f{btr)(8  =  ir,x) | 
vs  x.  The  quantity fibgr) ( ir,x )  is  as'in  Eq.  ( 1 )  but  with  A„ 
replaced  by  A  („bs'\  where  A  (rbgr),  for  most  metal  shells  in 
water,  is  the  rigid  one,  given  by:  A  f*''  =  —  j'„(x)/h  V  r(jc). 
We  have  often  denoted  these  plots  of  isolated  resonances  as 
the  “residual  responses”  since  these  residuals  constitute 
what  is  left  of  the  form  function  after  the  suitable  back¬ 
grounds  have  been  suppressed. 

(iii)  The  angular  scattering  pattern  at  a  fixed frequency: 
We  generate  and  display  |  fK  (6,  x  =  x0)|  vs<9.  Ifx0  isnota 
resonance  (i.e.,  a  real  root  of  the  real  characteristic  equation 
Ff'(x)  =  Re[zJ,u(x)  ]  ~  '),  then  the  resulting  plot  will  be 
complicated  and  uninformative.  Ifx0  is  a  (real)  resonance, 
then  the  associated  plot  displays  simplifying  and  informative 
features.*-9 

(iv)  The  (bistatic)  angular  scattering  pattern  at  a  fixed 
frequency,  after  background  subtraction:  The  pertinent 
quantity  here  is 

!/*  (<9.  x  =  x0)  ~r*Kr)(G,  x  =  x0)|  vs  0. 

The  quantity  f'b*r’(0.x() )  is  as  in  Eq.  (1),  but  with 
x  =  x0  =  fixed,  and  A„  replaced  by  A  \bgr\  as  in  (ii).  What 
was  done  above  for  the  summed  form  function  |  fx  (0,x)j 
can  also  be  done  for  the  partial  waves  contained  within  it 
[viz.,  \fn(0,x)\,  for  n  —  0,  1,  2, ...].  This  type  of  partial- 
wave  analysis  is  very  useful  and  informative. 

(v)  Generation  of  the  backscattered  partial  waves:  Here 
we  compute  and  display  |  f„(6  =  ~,x)  |  vsx.forn  =  0, 1,2, 
... .  Again,  Pn  (cos  sr)  —  (  —  1 )". 

(vi)  Resonance  isolation  within  each  mode  or  partial 
wave  of  index  n:  Here  we  compute  and  display 
|  f„(8  =  irpt)  —flnbt'>(8  =  irp c)|  vs x  in  the  backscattering 
direction  tr,  for  «  =  0, 1, 2, ... .  Hertf{b*')(TTrx)  is  as  in  (ii), 
with  oo  replaced  by  the  index  n. 

(vii)  The  (bistatic)  angular  pattern  of a  single  mode  n  at 
fixed  frequency:  The  pertinent  quantity  is  |/„(0,  x  —  x0 )  | 


vs  6,  for  each  value,  n  =  0,  1,2 . 

(viii)  The  angular  scattering  pattern  of  modal  reson¬ 
ances:  In  this  case,  we  generate  and  display 
I /„(&  *  =  x0)  —f''«*r)(0<  x --x0)|  vs  6,  for  n  =  0,  1,  2, 
... .  Here/,','',:')(7r,x0 )  is  as  in  (ii),  with  <x  replaced  by  the 
index  n.  If  the  fixed  x0  is  the  /  th  resonance  within  the  mode  n, 
(viz.,  x„(),  then  the  plot  will  be  proportional  to  the  angular 
pattern  of  a  Legendre  polynomial  |P„  (cos0)|,  and  will 
have  2 it  lobes  of  different  sizes.  In  the  cylindrical  case,  this 
plot  is  proportional  to  [cos  nd  [,  and  it  still  has  In  lobes  now 
of  equal  sizes.  Such  a  plot  resembles  a  “daisy”  or  a  “rosetta.” 
It  is  a  rhodonea  of  In  petals. 

It  is  often  instructive  to  generate  and  display  calcula¬ 
tions  versus  the  modal  index  n,  for  fixed  values  of  6  and/or  x. 
A  number  of  relevant  cases  are  as  follows. 

( ix)  The  response  surface:  Here  we  compute  the  isolated 
resonances  obtained  after  the  subtraction  of  the  appropriate 
background,  and  display  them  versus  n  and  x,  in  an  isomet¬ 
ric  view.  The  calculation  is  done  in  the  backscattering  direc¬ 
tion  0  =  sr  and  the  pertinent  display  is 
\fK  (6  =  tr,x)  —  fib*r' (0  —  tr,x)\,  versus  n  and  x.  This 
could  also  be  generated  in  other  directions. 10 

(x)  The  spectrogram  versus  mode  order  n,  at  fixed  fre¬ 
quencies:  Here  we  generate  and  display 

\fK  (0  =  tr,x  —  x0 )  -/(i*r)(m,x0 )  |  vs  //.  The  frequency  x0 
can  be  chosen  to  be  a  resonance  frequency  x„,  of  the  body,  or 
not.  This  is  a  slice  of  the  response  surface  in  item  (ix)  at 
x  =  x0.  Analogously,  the  spectrogram  in  item  (ii)  is  a  slice 
cf  the  response  surface  in  item  (ix)  at  each  integer  value  of  n. 

( xi )  The  isolated  medal  “residuals" displayed  versus  n  at 
fixed frequencies:  This  is  analogous  to  item  (x)  but  for  each 
partial  wave  or  mode  n.  We  generate  and  display 
I  fn(0~  ~,x0 )  ~fiber)(0 -  tt,x0 ) |  vs  n,  at  fixed  x0,  in  the 
backscattered  direction  0  —  sr.  Here,  x0  can  be  a  body  reso¬ 
nance,  xn„  or  not.  This  plot  shows  the  strength  with  which 
each  mode  n  contributes  to  the  resulting  residual,  at  each 
chosen  frequency  x0. 

(xii)  The  modal  backgrounds  in  mode-order  space  n: 
These  often  informative  plots  show  \f[bsr)(8  =  n\x0 )  |  vs  n, 
at  various  fixed  x  =  x0.  Clearly,  items  (xi)  and  (xii)  can  be 
used  to  displa-,  the  separation  of  “backgrounds”  and  “reson¬ 
ances”  in  the  mode-order  domain  n.  As  a  companion  to  item 
fvi),  j  f[bgr)(0  =  r,.x)|  can  also  be  displayed  versus  x  for 
n  =  0,  1,2,  ...  and  then  such  a  display,  together  with  item 
(vi),  exhibits  the  separation  of  “backgrounds”  and  “reson¬ 
ances”  in  the  frequency  domain,  which  is  standard  since  the 
eaily  days  of  the  development  of  the  RST.8"10 

We  have  always  distinguished  the  complex  eigenfre- 
quencies  x„„  roots  of  D„(x)  =0  in  the  complex  x  plane, 
from  the  real  resonances  x„f,  which  are  roots  of  leal  charac¬ 
teristic  equations,  which  in  this  case  are  of  the  type: 
Ff'(x)  =  Refz^’fx)]  “  ’.  It  isat  these  resonances x„,  that 
peaks  and/or  dips  appear  in  the  plots  of  the  form  functions 
versus  the  real  variable  x,  and  that  exactly  n  +  1/2  wave¬ 
lengths  fit  the  shell’s  circumference.  It  is  common  to  also 
consider  an  alternative  representation  in  the  complex  fre¬ 
quency  plane  x,  and  display  the  complex  ; -genfrequency 
there.  Both  representations  have  their  advantages.  The  com¬ 
plex  eigenfrequencies  of  a  metal  shell  in  water  tend  to  split 
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into  two  great  groups.  One  big  set,  clustered  just  beneath  the 
Re  x  axis,  is  associated  with  the  elastic  composition  of  the 
body.  The  other  large  set,  arranged  in  concentric  quasise¬ 
mielliptical  loci,  rnd  located  deeper  into  the  lower  half  of  the 
x  plane,  is  associated  with  the  body  shape.  To  obtain  these 
eigenfrequencies  (which  approximately  have  the  resonances 
?.s  real  parts)  one  only  needs  to  solve  an  equation  such  as 
Dn  (x)  =  0  foi  its  roots  xn„  and  it  is  obvious  that  this  can 
always  be  done  without  ever  having  to  subtract  “back¬ 
grounds”  as  one  must  do  if  we  look  at  the  frequency  depen¬ 
dence  of  the  form  function  after  the  subtraction  of  the  back¬ 
ground  [cf.  (ii)  ] .  The  point  here  is  that  once  the  resonances 
are  found  by  either  approach  (viz.,  either  subtracting  the 
appropriate  background,  or  without  ever  having  to  talk 
about  backgrounds  by  solving  for  the  roots  of  a  characteris¬ 
tic  equation),  they  can  be  used  to  construct  dispersion  plots 
for  the  phase  or  group  velocities  of  the  (surface)  elastic 
waves  that  the  resonances  generate  around  the  scatterer.  If 
we  obtain  the  roots  of  D„  (3c)  =  0  in  the  (complex)  x  p'.ine, 
then  their  real  parts  turn  out  to  be  close  to  the  resonances, 
and  their  imaginary  parts,  to  their  widths.  If  the  resonances 
are  isolated  by  means  of  the  spectrogram  in  item  (ii),  then 
their  widths  can  be  read  from  the  plot.  We  repeat  that  the 
resonances  thus  found,  exactly  coincide  with  the  real  roots 
x„,  of  the  (real)  characteristic  equation 
Ff'(x)  =  Re[z‘0(x)]  The  fact  that  tne  Rex,,,  are 
only  approximately  equal  to  thex„,  does  not  imply  arbitrari¬ 
ness  in  the  way  the  x„,  are  determined. 

(  xiii )  Dispersion  curves  for  the  phase  velocities  cf,  of  the 
surface  waves:  They  are  found  from  the  x„,  by  the  relation: 
cf/c,  =  x„,/(n  +  1/2). 

(xiv)  Dispersion  curves  for  the  group  velocities  c'f  of the 

surface  waves:  They  are  found  from: 

c{*Vc,  =  l/Re[cfo,(x)/dx],wherec,  is  the  sound-speed  in 
the  outer  medium.  Analogous  expressions  exist9  for  the 
phase  and  group  attenuation  “constants”  for  the  waves  re¬ 
volving  around  the  scatterer. 

(xv)  Results  in  vacuum  versus  fluid-loaded  results:  All 
the  results  discussed  above  are  usually  found  with  the  exact 
formulation  in  Refs.  5  and  6,  which  accounts  for  the  influ¬ 
ence  of  the  outer  and  inner  fluids.  A  formulation  for  shells  in 
a  vacuum  could  be  extracted  from  the  above  formulation 
simply  by  setting />,  =  0  andp?  =  0.  This  would  imply  that 
the  elements  du,  di6,  and  Af  (of  Ref.  6)  vanish.  Such  a 
result  would  be  unrelated  to  a  scattering  situation  and  would 
be  useful  only  to  describe  the  vibration  of  unloaded  elastic 
bodies.  In  this  case,  it  follows  that 

An(x)=A^'\x) 

=  fn(x)  _ _ 

hlf(x)x>l 

J/"e-tecosjc  (n  =  0,2,4,...) 

1  — i*e~feSinx  (rt=  1,3,5,...). 

These  shell  results  account  for  the  precise  curvature  of  the 
shell  in  a  natural  way,  whether  it  is  fluid  loaded  or  in  a  vacu¬ 
um.  It  is  obvious  that  ifsomeshell  feature  such  as  a  branch  of 
a  dispersion  plot  for  a  phase  velocity  were  to  be  displayed 
accounting  for,  and  then  ignoring,  the  fluid-loading  effects, 


the  resulting  curves  would  differ.  The  two  curves  would  then 
exhibit  a  discontinuous  “kink”  at  the  point  wherethe  fluid- ' 
loading  condition  changes. 

(xvi)  Zero-curvature  results:  There  have  been  attempts 
at  interpreting  features  in  the  scattering  cross  sections  of 
shells  by  means  of  results  derved  for  flat  plates.  Such  at¬ 
tempts  may  have  been  motivated  by  a  desire  to  see  how  a 
simpler  n.odel  could  explain  a  more  complex  situation  that  it 
was  not  originally  meant  to  explain.  It  is  also  conceivable 
that  they  may  have  been  due  to  the  difficulty  involved  in 
implementing  the  more  general  formulation.  At  h’gh-fre- 
quencies,  it  is  possible  to  estimate  the  phase  velocities  for  a 
shell  (fluid  loaded  or  in  vacuum)  by  their  values  obtained  for 
a  plate,  by  means  of  the  conversion: 


x  =  A-,0  — 


277 


(  fh)  n.n 

1  -b/a  ' 


(4) 


This  relation  converted  the  frequency-thickness  products  at 
which  certain  modes  (or  branches)  occurred  in  a  flat  plate, 
to  the  kta  values  for  which  “corresponding”  effects  were 
present  for  shells.  It  is  at  high  frequencies  that  the  shell  °ize 
becomes  large  compared  to  a  wavelength,  and  that  curvature 
effects  become  negligible.  The  approximation  in  Eq.  (4)  fails 
at  low  frequencies.  The  present  implementation  of  shell  for¬ 
mulations  (cf.,  Refs.  5  and  6),  is  valid  at  all  frequencies,  and 
thus,  it  brings  out  the  features  present  in  all  spectral  regions. 

In  what  follows,  we  will  always  use  our  shell  formula¬ 
tion  (Refs.  5  and  6).  It  will  always  account  for  exterior/in¬ 
terior  fluid-loading  effects,  as  well  as  for  the  precise  shell 
curvatures.  The  snell  deformations  are  described  by  the  (ex¬ 
act)  three-dimensional  equations  of  elasticity — as  opposed 
to  the  many  existing  “shell  theory”  approximations — and  it 
will  nowhere  need  estimates  based  on  flat  -plate  theories  or 
effects.  Hence,  our  results  will  be  valid  for  all  shell  thick¬ 
nesses,  curvatures,  and  in  all  frequency  bands,  whether  there  j 
is  fluid  loading  or  not. 

In  connection  to  the  fluid-loading  effects,  it  has  been 
repeatedly18  stated  that  the  dispersion  curve  for  the  zeroth- 
order  antisymmetric  Lamb  wave,  A0  propagating  in  a 
spherical  shell  in  a  fluid  starts  to  be  excited  for  Cp>c, .  For 
cp  <  c,  there  seems  to  be  no  A0  branch,  and  what  seems  to 
exist  then  is  a  section  of  the  dispersion  curve  corresponding 
to  the  “fluid-loaded”  case,  sometimes  called  the  A  branch. 
This  segment  seems  to  be  a  discontinuous  or  disjointed  ex¬ 
tension  of  the  A0  curve  toward  the  low-frequency  end  of  the 
spectrum.  See  below. 

In  connection  to  curvature  effects  it  also  follows  that  at 
low  frequencies  not  only  are  there  resonance  effects  due  to 
the  A0,  A,  and  S0  waves,  but  also  to  an  additional  curvature 
wave.  This  curvature  wave  is  missed  if  plate  theory  approxi¬ 
mations  are  used  to  predict  resonance  locations.  This  wave 
causes  the  breathing  mode  vibration  of  the  shell,  and  it  is  the 
shell  counterpart  of  the  interface  wave  responsible  for  the 
giant  (monopole,  n  =  0)  resonance91’  of  an  air  bubble  in  wa¬ 
ter.  For  a  shell,  this  resonance  is  not  giant  but  just  large,  and 
it  consists  of  about  half  a  dozen  multipole  modal  compo¬ 
nents  (as  can  be  seen  in  Fig.  3 )  or  its  enlargements.  (See  the 
discussion  of  Fig.  3,  below.)  Hence,  at  low  frequencies,  the 
general  exact  three-dimensional  shell  formulation  presented 
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in  Refs.  5  and  6,  and  implemented  here,  is  the  only  tool  that 
can  accurately  examine  arbitrary  fluid-loading,  curvature, 
and  resonance  effects  in  any  frequency  band  of  the  form 
functions  of  submerged  shells. 


II.  NUMERICAL  RESULTS  AND  DISCUSSION 

A  partial  analysis  of  figures  analogous  to  Figs.  1-3  has 
already  appeared5,6,18-20'23  in  various  earlier  works.  The 
complete  picture  that  emerges  from  in-depth  studies  of  these 
results  follows  in  very  abridged  form,  below.  With  the  excep¬ 
tion  of  our  own  work,5,6"8  the  form  functions  have  always 
been  analyzed  in  relative  narrow  bands  not  exceeding  A,  a 
values  of  100.  We  have  now  extended  these  bands  up  to 
A,  a  =  500,  which  allows  us  to  see  some  effects  never  before 
published,  or  even  noticed. 

Figure  1  shows  the  form  function  of  an  air-filled  steel 
shell  immersed  in  water.  The  relative  shell  thickness  is 
h/a  =  5%.  The  properties  of  steel  are:  cdi  =  5.95X105 
cm/s,  c3,  =  3.24  X  105  cm/s,  /?2  =  7.7  g/cm\  and  those  of 
water  art  c,  =  1.4825X  105  cm/s,/?,  =  1  g/cm3.  The  prop¬ 
erties  of  the  internal  air  are:  =  0.0012  g/cm3  and  c5 
=  0.344  X  105  cm/s.  Figure  2  displays  the  residual  response 
obtained  by  subtracting  the  rigid  background  over  the  entire 
spectrum,  for  the  same  shell  in  water  in  the  band 
0<A,a<500.  Figure  3  (top)  exhibits  the  form  function  of  a 
steel  shell  of  the  same  properties  as  before  but  of  thickness 
h/a  =  2.5%,  in  the  band  0<A,a<100.  The  bottom  plot 
shows  the  residual  response  resulting  when  we  subtract  the 
rigid  background  throughout  the  entire  spectrum.  Dips  in 
the  upper  plot  become  resonance  peaks  in  the  bottom  plot. 
These  curves  are  obtained  by  the  procedures  labeled  ( i )  and 
(ii)  in  Sec.  I,  and  the  (exact)  methodology  of  Refs.  5  and  6. 

Consider  Fig.  3  for  a  steel  shell  of  relative  thickness 
t~li  /a  =  2.5%.  There  is  a  clearly  visible  region  of  “strong 
flexures,”6  or  "bump,”  that  starts  to  develop  at  x~t  ~  *,  or 
in  this  case,  for  x—40.  The  empirical  condition  xt~  1  is  a 
rule  of  thumb  that  has  emerged  from  observation  of  many 
cases.6  '8  (See,  for  example,  Fig.  4  and  its  discussion  below, 
fora?  =  5%  shell,  or  Figs.  1-7  of  Ref.  6.)  It  has  its  origin  in 
the  standard  coincidence2'  condition  <?  =  c,  for  the  An 
wave.  The  frequency  a,  for  which  cp  —  c,  (see,  for  example, 
the  dispersion  plot  for  A0  in  Fig.  4).  seems  to  always  be  such 
that  xt—  1.  This  condition  depends,  in  general,  on  material 
parameters  through  a  proportionality  constant  C  (viz., 
Cxi~  1 ),  but  for  steel  Cs0.8,  and  thus,  it  still  follows  that 
xt  —  1 .  Thus,  as  a  rule,  the  region  of  strong  flexures  begins  to 
appear  at  the  coincidence  frequency.  In  addition  there  are 
resonance  peaks  spaced  at  regular  but  rather  broad  intervals 
throughout  the  whole  spectrum  except  around  the  coinci¬ 
dence518,19,21  frequency.  Their  separation  can  be  read  from 
Fig.  3  to  be  Ax  —  3.5  and  they  aredue  to  the wave,  which  is 
the  only  one  having  its  effects  felt  over  the  whole  spectrum.5 
Then,  near  coincidence  (viz.,  40<A,a<70),  the  spacing  of 
the  observable  resonances  becomes  smaller  (i.e.,  we  read 
Ax— 1.2)  and  nearly  uniform.  These  more  closely  spaced 
and  narrower  spikes  are  due  to  the  A0  wave,  which  is  almost 


FIG.  1.  Form  function  I/,  |  as  obtained  from  Eq  ( 1 )  for  a  spherical,  air- 
filled  steel  shell  of  relative  thickness /i  /a=\  —  b  /a  —  5%,  immersed  in  wa¬ 
ter  in  the  (nondimensional)  frequency  band  ()</;,  a<500. 


nondispersive,  has  lower  speed  (cf.  Fig.  4),  and  has  effects 
that  begin  to  be  felt  at  coincidence.5,6  We  also  note  some 
Iow-frequencv  resonance  features  at  A,  a  — 3,  which  are  due 
to  the  nonzero  curvature  of  the  shell.  For  the  case  of  very 
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FIG.  2.  Residual  (or  "resonance")  response  obtained  from  the  form  func¬ 
tion  in  Fig.  1  for  the  same  shell,  by  the  subtraction  of  the  (ordinary)  rigid 
background  over  the  entire  band  0 <  A,  o< 500. 
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Residual  Response 


Residual  Response  VS  Ka 


FIG.  3.  Form  function  (top)  and  residual  response  (bottom)  obtained  by  the  supression  of  the  rigid  background  for  a  steel  spherical  shell  of  thickness 
ft /as  2.5%  in  water.  These  plots  show  the  region  of  strong  flexures  around  the  coincidence  frequency,  for  which  cf  —  c, . 
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FIG.  4  Exact  dispersion  plots  for  the  phase  velocities  cp  of  the  first  few 
symmetric  and/or  antisymmetric  Lamb  waves  circumnavigating  a  steel 
spherical  shell  of  thickness  h  /a  —  5%,  in  water.  The  bottom  plot  is  an  en¬ 
largement  of  the  top  one  for  the  !ow-frequency  region  enclosed  in  the 
dashed-line  rectangle.  The  top  graph  is  the  exact  spherical  counterpart  of 
the  often  quoted  Fig.  26  of  Ref.  1 2.  for  plates. 


thin  shells,  these  resonance  features  merge  to  produce  the 
(monopoie)  breathing  mode  of  an  air  bubble  in  water,  which 
becomes  stronger  and  eventually  giant,  the  thinner  the  shell 
becomes.  Figures  1  and  2  also  reveal  a  high-frequency  fea¬ 
ture  (viz.,  A,  a  —  240)  that  had  gone  unnoticed  in  the  past  for 
shells  of  this  thickness  since  it  appears  for  k,a>  100 — a  re¬ 
gion  not  investigated  much  in  the  past.  For  thinner  shells, 
such  a  broad  resonance  occurs  at  even  higher  frequencies 
(viz.,  at  A,a~480  for  t  =  2.5%;  cf.  Fig.  9),  as  we  will  see 
later. 

Once  we  have  obtained  the  resonance  locations,  we  can 
use  them  to  construct  dispersion  plots  for  the  phase  veloc¬ 
ities  of  the  propagating  (Lamb)  waves.  Resonance  locations 
can  be  read  from  Fig.  2  or  3  (bottom),  or  even  from  Figs.  1 
and  3  (top),  since  in  these  instances  the  resonances  are  quite 
noticeable  even  before  isolating  them  by  background  sub¬ 
traction.  As  we  mentioned  in  Sec.  I,  item  (xii),  resonance 
locations  can  also  be  obtained  by  solving  for  the  (real) 
roots  of  characteristic  equations  such  as  F~  ‘(x) 
=  Refz^'fx:)]  - ',  and  then  substituting  them  into  the 
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well-known  relation  given  in  item  (xiii)  of  Section  I.  This 
alternative  procedure  does  not  require  background  subtrac¬ 
tion,  but  is  more  computationally  intensive.  In  either  case, 
the  result  of  the  analysis  for  a  /  =  5%  steel  shell  in  water  in 
the  band,  0 ;u<500,  is  shown  in  Fig.  4  (top).  A  detail  of 
this  plot  in  the  low-freque.ic>  region  (da  bed  rectangle)  is 
shown  in  Fig.  4  (.bottom).  This  figure  contains  much  infor¬ 
mation  as  we  wiii  see  below.  As  memioned  in  item  (xvi), 
there  have  been  attempts  at  interpreting  the  form-function 
features  of  shells  in  terms  of  "corresponding”  results  for  flat 
plates.  Consider  the  often-quoted  rig.  26  of  Rrekhovskikh’s 
book  (Ref.  12)  which  displays  the  dispersion  plots  for  the 
phase  velocities  of  symmetric  Sj  and  antisymmetric  Aj 
(j  =  0,1,2...)  Lamb  waves  on  an  aluminum  plate  in  vacuum. 
Brekhovskikh  attributes  this  figure  to  Schoch25  who  in  turn 
attrib  tes  it  to  Firestone.26  We  have  recalculated  it  for  a 
tungsten  carbide  plate  (Fig.  4,  Ref.  5).  Many  authors  have 
generated  it  for  various  metals.' 4  That  figure  shows 
that  at  low  frequencies,  the  A0  branch  tends  to  zero  and  the 
S0  curve  approaches  a  small  constant  value  slightly  below 
the  value  cdi  /cs% .  Figure  4  is  the  counterpart  of  that  old  plot, 
now  for  air-filled  steel  spherical  shells  in  water.  In  the  band 
0  A,  a  <  500,  five  branches  enter  the  picture,  namely  A0 ,  A , , 
and  Cn,  S’,,  and  S\.  For  shells,  both  the  A0  and  the  50 
branches  attain  very  high  c  values  at  the  low-frequency  end 
ofthe  plot.  Having  Fig.  4,  it  is  no  longer  necessary  to  use  flat- 
plate  results  (such  as  Fig.  26  of  Ref.  12)  to  physically  inter¬ 
pret  resonance  features  in  the  backscattering  cross  sections 
of  fluid-loaded  elastic  shells.  In  fact,  Fig.  4,  which  corre¬ 
sponds  to  the  exact  solution  for  a  shell  treated  by  three-di¬ 
mensional  elastodynamics,  differs  substantially  at  low  fre¬ 
quencies  from  the  equivalent  (standard)  results 
obtained1213  from  plate  theories.  According  to  the  relation 
in  Sec.  I  (xiii),  the  phase  velocities  c7(x)  are  proportional  to 
the  resonances  xnl,  found  as  roots  cf  the  pertinent  character¬ 
istic  equation.  This  whole  procedure  is  entirely  contained  in 
the  frequency  domain  x,  and  is  urn  elated  to  approaches  de¬ 
veloped  for  the  mode-order  domain.  The  numbers  marked 
along  the  various  branches  ol  the  dispersion  curves  corre¬ 
spond  to  values  of  the  index  n  at  the  various  frequencies  at 
which  they  occur.  These  numbers  are  obtained  by  carrying 
on  a  partial-wave  analysis  of  tne  form  function  as  indicated 
in  items  (v)  and  (vi)  ofSec  I,  for  various  values  ofthe  mode 
index  n.  At  high  frequencies,  the  A0  branch  tends  to  the 
Rayleigh  speed  cR.  Figure  4  (bottom)  shows  that  the  A0 
branch  crosses  the  c  =  c,  level  at  the  coincidence  frequency, 
near  x  =  t  ~  1  ~20.  For  frequencies  below  coincidence,  the 
A0  branch  is  to  be  replaced  by  the  fluid-loaded  branch  A, 
which  here  is  determined  from  our  (exact,  3-D)  shell  formu¬ 
lation, 5,6  not  from  any  flat-plate  result.25,26  Only  in  this  fash¬ 
ion  can  the  shell  curvature  (as  well  as  fluid  loading)  effects 
be  preserved  and  studied.  The  (sometimes  called  fluid  load¬ 
ed  )  A  branch  is  not  shown  here  ( in  Fig.  4 ),  but  it  is  known  to 
remain  always  "subsonic”  c  <  c, .  The  term  "subsonic,”  or 
“supersonic”  (c>  c, )  often  used  to  describe  these  phase  ve¬ 
locities,  is  to  be  used  with  the  understanding  that  there  are  no 
perceivable  transonic  effects  or  "sonic  booms"  occurring  as 
these  phase  velocities  exceed  the  sound  speed  in  the  medium. 
It  is  known  that  phase  v  :locities  in  analogous  electromagnet- 
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FIG  5.  Exact  dispersion  plots  for  the  phase  velocities  cr  of  the  Lamb  waves 
on  a  water-loaded  steel  spherical  shell  of  thickness  h  /a  =  2%  in  the  band 
0<A',tf<500.  Three  branches  of  the  dispersion  curves  (viz.,  those  corre¬ 
sponding  to  the  A0,  S0.  and  waves)  are  present  in  this  band.  The 
numbers  along  the  curves  correspond  to  values  of  the  index  n  at  the  various 
frequencies. 


ic  instances  can  even  exceed  the  speed  of  light  in  the  medium 
in  question.  The  S0  branch  is  the  exact  3-D  elasticity  coun¬ 
terpart  of  lunger’s27  upper  branch  of  eigenfrequencies  found 
by  means  of  Donnell’s  shell-theory  equations.28  At  low  fre¬ 
quencies,  the  Donnell  theory  results  do  not  differ  much  from 
the  exact  3-D  elasticity  results  given  in  Fig.  4,  for  that 
branch.  Note  that  the  exact  S0  branch  displayed  in  Fig.  4  has 
a  relatively  flat  shape  in  the  enlargement  (bottom),  but  as  it 
is  seen  in  the  top  graph  for  0</r,o<500,  at  the  higher  fre¬ 
quencies,  it  bends  down  and  asymptotically  approaches  cR , 
just  as  the  A0  branch. 

Figure  5  displays  analogous  results  for  a  spherical  steel 
shell  in  water  of  thickness  t  —  2%  in  the  band  0<A',  c<500. 
For  this  thickness,  only  three  branches  of  the  dispersion 
curves  (viz.,^0,  A,,  and  S0)  are  present.  The  others  appear 
at  still  higher  frequencies.  Figure  6  displays  analogous  re¬ 
sults  for  the  shell  in  the  same  band  when  the  thickness  is  now 
t=l%.  Only  the^0  and  S0  branches  are  present  in  this  case. 
All  that  was  said  for  the  thicker  shell  still  applies  for  Figs.  5 
and  6.  The  bottom  part  of  Fig.  6  shows  a  section  of  the  resid¬ 
ual  response  in  the  high-frequency  band  400  <A-,  a  <500  for  a 
r  =  1%  steel  shell  in  water.  In  general,  suppression  of  the 
rigid  background  does  not  exactly  isolate  the  resonances, 
everywhere,  in  a  completely  clean  and  clear  fashion .  The  last 
segment  of  Fig.  2  (400<£,c<500)  for  the  t  =  5%  shell 
shows  the  point  well.  Even  though  at  high  frequencies  the 
rigid  background  should7  become  “better,”  the  last  segment 
of  Fig.  2  shows  that  such  is  not  always  the  case.  Figure  6 
(bottom)  is  displayed  here  merely  to  show  that  if  a  "modi¬ 
fied”  background  is  used,  the  resonance-isolation  process 
that  results  is  very  clean  and  accurate.  The  modified  back¬ 
ground  responsible  for  Fig.  6  does  not  assume  the  scatterer 
rigid  as  in  many  earlier  RST  papers,  but  it  accounts  for  the 
matching  of  elastic  properties  of  the  scatterer  and  the  fluid 
loadings  taking  place  at  the  shell’s  surfaces.  Its  mathemat- 
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FIG.  6.  (i)  Exact  dispersion  plots ofthe  phase  velocities^  ofthe,40  and  S’,, 
Lamb  waves  (top)  propagating  in  a  thin  steel  spherical  shell  of  h  /a  =  1  %  in 
water  in  the  band  0<A1o<500.  These  are  the  only  two  branches  present 
within  this  band,  (ii)  Portion  of  the  residual  response  (bottom)  of  the  same 
shell  in  the  band  400<A-,aO00.  Notice  how  “clearly”  are  the  resonance 
peaks  isolated  in  this  plot  which  resulted  from  the  s'!1” -action  of  a  "modi¬ 
fied"  background. 


ical  derivation  is  lengthy  and  will  appear  elsewhere.  Simple 
observation  of  Fig.  6  is  the  best  evidence  of  its  appropriate¬ 
ness. 

Figure  7  shows  the  dispersion  plot  for  the  phase-velocity 
of  the  Aq  branch  for  a  steel  spherical  shell  of  t  =  1%  in  wa¬ 
ter.  It  coincides  almost  exactly  with  the^0  branch  shown  in 
Fig.  6  (top)  which  was  computed  using  the  exact  3-D  for¬ 
mulation  of  Refs.  5  and  6.  However,  it  was  calculated  hereby 
means  of  the  Donnell  shell-theory  equations38  used  by 
Junger  and  others.37  This  plot  is  simpler  and  faster  to  gener¬ 
ate  than  the  one  in  Fig.  6  (top). 

Figure  8  shows  the  form  equation  (top)  and  the  residual 
response  (bottom)  of  a  steel  shell  of  t  =  1%  in  water.  The 
curves  are  found  in  the  band  300<&,fl<500,  and  the  lower 
ones  have  made  use  of  the  "modified”  background  men¬ 
tioned  above.  While  the  form  function  (top)  exhibits  some 
peaks  and  some  dips,  the  “residual”  response  displays  clear 
resonances  everywhere.  In  fact,  as  stated  above,  these  residu¬ 
als  are  the  isolated  resonances,  as  predicted  by  the  present 
approach. 

Figure  9  is  analogous  to  Fig.  8,  but  the  steel  shell  now 
has  thickness  t  =  2.5%.  The  novelty  is  the  appearance  of  a 
large  broad  peak  near  A:,o~486.  In  fact,  it  extends  from 
kj  a  =  470  to  5 10.  This  feature  was  seen  earlier  at  k ,  a~240 
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FIG.  7.  Detail  of  the  dispersion  plot  for  the  phase  velocity  cf  of  the/f0  Lamb  wave  on  a  1%  steel  shell  in  water,  in  the  band  0<kla<  100. 


FIG.  9.  Same  as  Fig.  8  but  now  for  a  steel  shell  of  thickness  2.5%.  The  main 
feature  observable  in  this  case  is  a  large  and  wide  resonance  peak  near 
FIG.  8.  Form  function  ( top)  and  residual  response  (bottom)  for  a  1  %  steel  k,  a~  480.  This  is  caused  by  a  thickness  resonance  effect  resulting  when  the 

shell  in  water  in  the  band  300<*,  o<500.  We  underline  the  clean  and  clear  frequency-thickness  product  for  the  shell  equals  either  the  phase  velocity  of 

nature  of  the  isolated  resonances  in  the  bottom  plot,  obtained  via  the  same  the  pertinent  Lamb  wave  propagating  along  the  shell,  or  half  the  shell’s 

“modified"  background  used  in  Fig.  6  (ii).  dilatational  wave  speed. 
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FIG.  10.  Phase  velocities  c,  of  the  three  Lamb  waves  (viz.,  A0%  50,and/4,) 
that  exist  for  kt  a< 500  for  a  steel  shell  of  thickness  h  /a  =  2.5%  in  water, 
displa>ed  in  the  band  300<.kla<S00.  For  An,  at  a — 480,  we  read 
cfi  ~2.87  Km/s. 


for  a  5%  shell  in  Figs.  1  and  2.  Figure  10  displays  dispersion 
plots  for  the  phase  velocities  of  the  first  three  Lamb  waves 
A0,  S0,  and/4,  in  the  band:  300<A,  a<500.  We  note  that  near 
480  the  cp  of  the  A0  branch  is  —2.87  Km/s.  This  new  large 
peak  is  caused  by  a  thickness  resonance  effect  present  in  the 
shell  which  we  analyze  below.  This  peak  could  be  caused  by 
the  fact  that  the  phase  velocity  of  the  propagating  wave 
equals  the  frequency-thicknesyjroduct  (viz.,  for  cp  =fh  ). 
This  condition  is  not  due  to  the  coincidence  phenomenon, 
although  it  resembles  it  in  form.  In  fact,  this  resonance  con¬ 
dition  could  be  simply  stated  as:  h  =  A,  where  A  =  cp/f  is 
the  wavelength  of  the  A0  wave  in  the  shell,  of  phase  velocity 
cp.  Shell  motions  ir.  this  situation  are  of  a  bending  or  flexural 
nature,  and  follow  a  direction  through  the  thickness,  normal 
to  the  fluid-solid  interface,  as  one  would  observe  at  the  inci¬ 
dence  point  on  the  shell.  Such  condition  serves  to  predict  the 
spectral  location  of  the  large  thickness  resonance,  since 


x  =  k^a- 


lira  j._  lira  cp  _  2tr  cp 
c,  h  c,  h/a 


(5a) 


*- 1  '-i  «  *■  i 

Substitution  of  the  numerical  values  yields 
2x3. 14X2. 87X  10s 


1.4825  X105  X  0.025 


=480, 


just  as  seen  in  Fig.  9.  There  is  a  “principle  of  transverse 
resonance”1  U2,14,24  that  could  also  account  for  this  peak, 
and  which  is  usually  written  as  h  =  mA,/2  for  m  =  2, 4,  6, 
....  where  A,  is  the  shear  wavelength  in  the  layer.  For  m  —  2, 
this  condition  is  A  =A,;  (or  equivalently,  k,h  —  mir),  which 
differs  from  the  previous  resonance  condition  in  that  we  now 
have  A,  rather  than  A  in  the  right  side.  This  leads  toe,  =Jh 
and  to  an  equation  to  predict  the  peak  location  which  is 


x  —  ■ 


s 


h/a  c, 


(5b)  , 


which  differs  from  Eq.  (5a)  in  that  cp  is  now  replaced  by  c,. 
There  is  also  an  analogous  “principle  of  longitudinal  reso¬ 
nance”1  U2J4,24  usually  written  as  h  =  n{Ad/2)  for  n  —  1,2, 
3,... ,  where  Ad  is  thedilatational  or  longitudinal  wavelength 
in  the  layer,  which  is  easily  seen  to  be  equivalent  to 
kdh  =  nir.  For  n  =  1,  it  leads  to  cd/2  —  jh,  or  simply, 
h  —  Ad/2.  The  peak  location  could  be  now  predicted  by 


2  v  cd 
h/a  2c, 


(5c) 


which  differs  from  Eq.  (5a)  in  that  cp  has  been  replaced  by 
cd/2.  Substituting  the  values  of  the  material  parameters  into 
Eqs.  (5b)  and  (5c)  yields 

6.28x3.24X10'  .,a 

x  = - -  -  549, 

0.025  X  1.4825  X  105 

for  Eq.  (5b),  and  x  =  504  for  Eq.  (5c).  Figure  9,  which 
displays  our  numerical  prediction,  shows  the  broad  thick¬ 
ness  pea:-  extending  within:  470<x<  5 10,  with  the  top  of  the 
peak  at  —480.  This  rules  out  the  prediction  of  Eq.  (5b), 
which  is  — 13%  too  large.  However,  it  implies  that  both  Eqs. 
(5a)  and  (5c)  predict  a  location  for  the  thickness  resonance 
peak  which  is  in  agreement  with  our  calculations  in  Fig.  9  for 
a  steel  shell.  In  the  absence  of  other  evidence,  it  seems  that 
the  prediction  of  the  main  peak  displayed  in  Fig.  9  is  best 
given  by  Eq.  (5a);  however,  that  of  Eq.  (5c)  4s  only 
1 8/504  £g  3.5%  away  from  the  other,  and  thus,  the  differ¬ 
ences  are  small.  They  are  close  enough  to  each  other  to  be 
both  “confirmed”  by  the  calculations  we  have  performed  for 
several  other  thicknesses  and  materials  which  will  be  shown 
elsewhere.  At  present,  it  is  not  clear  if  the  location  of  this 
large  thickness  resonance  peak  should  be  determined  by  that 
of  its  central  lobe  of  highest  amplitude,  or  by  that  of  one  of 
the  particular  subpeaks  that  form  it  (see  Fig.  9).  This  point 
will  be  further  studied  by  us  in  the  future.  For  the  case  in  Fig. 
2,  the  phase  velocity  is  slightly  less  than  the  values  found 
here,  so  the  location  of  the  peak  in  Fig.  2  is  about 
x=kf  o=236,  as  seen  there.  Finally,  the  resonance  fora  1% 
steel  shell  in  water  is  enlarged  in  Fig.  11  in  the  band 
90<  A,  a  <  1 50.  It  is  in  this  band  that  the  strong  flexures  begin 
to  appear  at  the  coincidence  frequency  x=  100.  This  figure 
was  already  examined  in  Fig.  3(c)  of  Ref.  6;  however,  now 
the  various  peaks  in  Fig.  1 1  have  been  identified  with  the 
various  mode  orders  n  that  cause  them.  This  is  done  by  a 
partial-wave  analysis  as  described  in  items  (v)  and  (vi)  of 
Sec.  I,  and  used  in  Figs.  4-6.  This  partial-wave  analysis 
shows  each  one  of  the  individual  contributions  to  the  ensem¬ 
ble  of  peaks  appearing  in  Fig.  1 1,  and  permits  one  to  identify 
them  in  that  fashion.  Figure  1 1  displays  the  set  of  resonances 
numbered  n  =  107  to  126,  which  are  separated  a  narrow 
distance  A*—  1.2,  and  are  due  to  the/4-  .,  ..hich  in  the 
present  exact  formulation,  always  acounts  for  fluid  loading 
and  shell  curvature.  It  also  displays  another  resonance  set, 
spaced  a  wider  amount,  Ax— 3.5,  toward  the  right  end  of  the 
plot,  which  is  due  to  the  S0  wave,  originally  studied  by 
Junger27  by  means  of  a  shell-theory  approximation  analo¬ 
gous  to  the  one  we  used  in  Fig.  7.  Coincidence  occurs  near 
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k,  a~  108.  For  a  t  =  1%  shell,  it  should  appear  at  x  =  k ,  a 
=  t  ~ 1  =  100  but  as  we  mentioned  above,  this  relation  is  to 
be  modified  by  a  multiplicative  factor  that  depends  on  mate¬ 
rial  parameters  and  that  is  close  to  unity.  The  left  portion  of 
Fig.  11,  up  to  ka~  125  actually  shows  the  superposition  of 
two  resonance  families.  The  flexural  (antisymmetric)  reso¬ 
nance  family  begins  to  appear  slightly  below  the  coincidence 
frequency,  at  which  point,  a  (forthcoming)  partial-wave 
analysis  of  the  figure  shows  that  the  individual  resonance 
contributions  are  quite  broad  when  displayed  versus  kc. 
These  flexural  contributions  are  coherent  and  add  in  phase 
below  coincidence,  but  are  out  of  phase  above  it  so  that  an 
envelope  type  of  function  becomes  noticeable.  This  envelope 
makes  up  the  “bump”  characteristic  of  this  neighborhood  of 
the  coincidence  frequency.  Superimposed  on  this  envelope  is 
a  narrow  resonance  family  due  to  a  surface  wave  analogous 
to  the  Stoneley  wave  at  the  interface  of  an  elastic  plate  fluid- 
loaded  on  only  one  side.  Alternative  names  for  these  reson¬ 
ances  are  “leaky  pseudo-Stoneley  resonances,”  and  also 
“Junger-type  resonances,”  and  they  appear  near  coinci¬ 
dence.  At  their  inception,  and  for  the  steel  in  Fig.  11,  they 
have  a  phase- velocity  of  0.88c, ,  which  ultimately  increases 
to  c,  at  coincidence,  at  which  point  they  are  essentially  dissi¬ 
pated.  In  the  vicinity  of  the  coincidence  frequency  where 
they  occur,  these  two  families  seem  to  be  the  same,  but  they 
really  have  different  origins.22  This  point  can  be  verified  by 
calculating  the  responses  of  shells  to  incident  pulses  at  coin¬ 
cidence,  in  the  late-time  region.  This  excludes  the  broader 
flexural  family  (A0 )  which  dissipates  faster  in  time,  and  per¬ 
mits  the  extrac.ion  of  the  associated  group  velocities  from 
the  resonance  half-widths,  of  the  Junger-type  family,  as  was 
reported  earlier.29  A  more  detailed  description  and  study  of 
the  partial-wave  analysis  that  leads  to  Fig.  1 1  and  these  re¬ 
sults,  will  appear  elsewhere.  Items  (viii)-(xii)  of  Sec.  I,  al¬ 
though  not  used  here,  have  been  used  in  past  resonance  stud¬ 
ies®"10,  and  will  be  used  in  future  ones. 


III.  CONCLUSIONS 

The  present  paper  examines  the  features  in  the  scatter¬ 
ing  cross  sections  of  air-filled  spherical  steel  shells  in  water  in 
the  band  0<A',  c<500.  The  resonance  features  are  then  asso¬ 
ciated  with  the  Lamb  waves  propagating  in  the  shell  materi¬ 
al,  and  with  other  external  circumferential  waves  propagat¬ 
ing  in  the  outer  fluid.  The  presence  of  circumferential  waves 
in  the  inner  fluid  does  not  manifest  itself  appreciably  in  the 
scattering  cross  sections.  Having  these  resonance  responses 
it  is  not  hard  to  work  backwards,9  and  be  able  to  identify 
various  properties  of  the  scatterer  such  as  its  overall  dimen¬ 
sion,  thickness,  and  material  composition.  For  example, 
knowledge  of  the  location  of  the  large  lobe  that  develops  at 
the  (nondimensional)  coincidence  frequency,  x,  for  which 
jcr~  1,  immediately  gives  the  relative  shell  thickness  t. 

Our  analysis  is  based  on  an  exact  three-dimensional 
elasticity  methodology  that  we  described  in  Refs.  5  and  6, 
initially  implemented  there,  and  continued  here.  This  meth¬ 
odology  is  ideal  to  reliably  analyze  any  scattering  problem 

J.  Acoust.  Soc.  Am.,  Vol.  89,  No.  4,  Pt.  1 ,  April  1991 


FIG.  1 1.  Detail  of  the  resonance  response  of  a  steel  spherical  shell  of  ( rela¬ 
tive)  thickness  h/a-~  1%,  air  filled,  and  submerged  in  water,  in  the  band 
90<£,a<150  in  the  vicinity  of  the  coincidence  frequency,  at  which  the 
strong  flexures  are  seen  to  develop.  The  resonance  peaks  are  labeled/num¬ 
bered  by  the  partial-wave  analysis  described  in  the  text 


pertaining  to  insonified  shells  fluid  loaded  on  their  inner  and 
outer  surfaces,  particularly  at  low  frequencies.  This  ap¬ 
proach  never  misses  any  of  the  many  effects  present  in  the 
scattering  process  whether  they  are  due  to  curvature,  fluid 
loading,  thickness  resonances,  or  to  choice  of  frequency 
band. 

There  is  no  need  to  analyze  or  interpret  resonance  fea¬ 
tures  in  shell  cross  sections  by  means  of  plate  waves  or  plate 
vibrational  theories,  whether  these  are  fluid  loaded  on  one  or 
two  sides,  or  placed  in  a  vacuum.  The  tendency  to  do  that  has 
been  quite  popular,  and  it  has  included  ourselves  (cf.  Fig.  4, 
Ref.  5 ) .  However,  at  present,  we  realize  that  it  only  obscures 
the  issues  and  introduces  errors.  It  is  clear  from  the  present 
study  that  only  the  3-D  elasticity  formulation  ol  liefs.  5  and 
6  accounts  for  all  possible  shell  effects.  Next  best,  is  the  use  of 
a  shell  theory  approximation27-28  which  will  be  ccurate  at 
least  up  to  the  coincidence  frequency.  As  we  saw  in  Fig.  7, 
such  an  approach  can  still  extract  the  curvature  effects  re¬ 
sponsible  for  the  high  values  of  the  phase  velocities  along  the 
Aq  and  S0  branches  at  low  frequencies.  (We  repeat  that  the 
branch  we  have  called  A0  already  contains  the  fluid-loading 
effects  within  it.)  The  basic  explanation  is  that  above  coinci¬ 
dence  there  is  an  A0  (Lamb)  wave  in  the  shell  and  also  a 
corresponding/^  wave  in  the  fluid,  which  is  transmitted  to  V 
by  the  motion  of  the  shell.  Below  coincidence,  the  A0 
( Lamb )  wave  in  the  shell  is  turned  off  or  ceases  to  be  excited, 
while  its  counterpart  in  the  fluid  does  not.  This  counterpart 
has  to  account  for  the  presence  of  the  surrounding  fluid, 
since  without  it,  it  would  have  no  place  to  propagate.  Thus 
what  has  been  termed  as  a  {plate)  A  wave  (that  “accounts” 
for  the  fluid  loading)  emerges  in  a  natural  way  from  the 
present  methodology  which  obviously  always  accounts  for 
the  fluid  loading  on  the  (curved)  shell.  We  remark  that  the 
present  exact  three-dimensional  methodology  is  also  avail¬ 
able  for  layered  shells.3  This  formul  .ion  also  leads  to  the 
generation  of  exact  dispersion  plots  for  the  phase  velocities 
of  the  various  types  of  Lamb  waves  propagating  on  fluid- 
loaded  spherical  shells.  This  generalizes  early  plate  results, 
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such  as  those  in  Fig.  26  of  Ref.  12  (or  Fig.  4,  Ref.  5),  and 
discourages  possible  future  attempts  to  interpret  resonance 
features  in  shell  cross  sections  by  plate  waves  or  results.  We 
have  favorably  compared  dispersion  plots  obtained  by  the 
present  exact  methodology  anu  those  based  on  the  Don¬ 
nell27  shell  equations  (Fig.  7).  We  have  studied  a  large-am¬ 
plitude  thickness  resonance  effect  present  in  shells  at  rela¬ 
tively  high  frequencies.  Since  calculations  are  rarely  carried 
out  to  very  high  frequencies  (viz.,  x~500  and  beyond),  we 
suspect  that  it  was  precisely  for  that  reason  that  it  had  re¬ 
mained  unnoticed.  The  present  broad  band  calculations  ex¬ 
hibited  here  have  brought  it  into  the  open.  The  location  of 
this  thickness  resonance  can  be  predicted  by  Eqs.  (5a)  and 
(5c),  which  rest  on  resonance  conditions  of  the  simple  type: 
h=A,orh=  hd/2,  respectively.  A  third  possible  resonance 
condition  (viz.,  =  h  )  examined  in  this  regard,  does  not 
seem  to  agree  with  the  nume~cal  calculations  we  have  per¬ 
formed  and  displayed  in  Fig.  ).  However,  these  calculations 
(cf.  Fig.  9)  do  confirm  the  predictions  of  both  Eqs.  (5a)  and 
(5c)  within  acceptable  accuracy.  It  also  follows  that  if  calcu¬ 
lations  were  to  be  performed  in  even  broader  frequency 
bands  (i.e.,  x>  500)  than  shown  here,  one  would  be  able  to 
note  higher-order  thickness  resonance  effects  of  the  types 
discussed  here,  caused  by  the  overtones  associated  with  mul¬ 
tiples  of  the  appropriate  wavelengths.  We  have  illustrated 
(in  Fig.  6,  bottom)  that  the  resonance-isolation  process  typi¬ 
cal  of  the  RST  can  be  very  clearly  and  accurately  accom¬ 
plished  over  the  entire  spectrum  by  the  suppression  of  a 
“modified”  background  that  will  be  the  subject  of  a  later 
detailed  study.  Furthermore,  by  means  of  a  partial-wave 
analysis,  we  have  identified  the  basic  types  of  resonance  fea¬ 
tures  present  in  the  resonance  response  of  a  submerged  steel 
shell  in  the  vicinity  of  its  coincidence  frequency,  and  estab¬ 
lished  their  connection  with  the  underlying  surface  waves. 
We  remark  in  closing  that  the  features  in  the  cross-section  of 
any  submerged  shell,  as  obtained  from  an  incident  cw  wave¬ 
form,  are  related  to  those  occurring  when  the  incidence  is 
pulsed.  For  example,  an  energetic  sinusoidal  pulse  of  short 
duration  incident  on  a  shell  will  backscatter  a  pressure  re¬ 
sponse  that  will  reproduce  the  (steady  state)  form  function 
I  /«  (*)  I  that  was  obtained  for  cw  incidences.  For  this  to 
occur,  the  carrier  frequency  x0  of  the  incident  pulse  should 
coincide  with  one  of  the  shell  resonances  x{r).  This  replica¬ 
tion  of  the  response  will  be  more  evident  for  frequencies  be¬ 
low  that  carrier  frequency.  Conversely,  for  long  incident 
pulses  (say,  having  a  delta  function  centered  around  x0  for 
spectrum),  the  amplitude  of  the  scattered  pressure  field  will 
reproduce  the  portion  of  |  /K  (x)|  near  a  shell  resonance 
provided  that  x0  =  x{r).  In  other  words,  long  incident  cw 
pulses  with  x0  =  x<r),  will  excite  individual  resonances  x(r> 
in  the  target.  These  ideas  can  all  be  extracted  from  the  filter- 
type  relation 

rp^)=~  f  +  °7«  (x')Pinc(x')e-^dx‘,  (6) 

2tt  J- « 

which  gives  the  scattered  pressure  in  the  pulsed  case,  in 
terms  of  the  form  function  in  the  steady-state  case,  and  the 
spectrum  of  the  incident  pulse.  In  this  later  case,  knowledge 
of  the  transient  vibration  of  the  shell30  becomes  critical. 
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